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f**^ _ Abstract. We investigate the particle trajectories in solitary waves with vorticity, where the 

, vorticity is assumed to be negative and decrease with depth. We show that the individual particle 

, ^ ■ moves in a similar way as that in the irrotational case if the underlying laminar flow is favorable, 

' that is, the flow is moving in the same direction as the wave propagation throughout the fluid, and 

show that if the underlying current is not favorable, some particles in a sufliicicntly small solitary 
wave move to the opposite direction of wave propagation along a path with a single loop or hump. 
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1. Introduction 

The water-wave problem concerns the gravity-driven flow of a perfect fluid of unit 
density; the flow is bounded below by a rigid horizontal bottom {Y = —d} and above by a 
free surface {Y = r]{X,t)}, where rj depends upon the horizontal spatial coordinate X and 
time t. Steady waves are waves which propagate from left to right with constant speed c 
and without change of shape, so that ri(X,t) = rj{X — ct). Solitary waves are steady waves 
which have the property that r](X — ct) as X — ct ^ ±00. We consider in this paper the 



. particle trajectories in a fluid as a solitary wave propagates on the free surface, assuming 



that the flow admits a negative vorticity decreasing with depth. 

There have been a series of works concentrating on the study of solitary waves, in 
the setting of both irrotational flows [H [H O [6l [7] etc., and rotational flows which become 
active only in the last few years. One of the interests in the above works is the description 
of individual particle path in the fluid. In irrotational flow, particle paths underneath a 
solitary wave are investigated in [U [6] , both for the smooth solitary wave and the solitary 
^ . wave of greatest height for which the crest is a stagnation point. It was shown in [6j that 

H ! in an irrotational solitary water wave, each particle is transported in the wave direction 

but slower than the wave speed; as the solitary wave propagates, all particles located ahead 
of the wave crest are lifted while those behind have a downward motion, with the particle 
trajectory having asymptotically the same height above the flat bed. For rotational flow, 
some results on particle paths under periodic waves are obtained (cf. |1H \TU[ \W[ [T8] for 
instance), following a series of works on the corresponding results in irrotational case (see 
[a El [13] etc). 

Recently, rigorous existence results on small-amplitude solitary waves with arbitrary 
vorticity distribution were obtained in [TJ] and in [12], using generalized implicit function 
theorem of Nash-Moser type and spatial dynamics method, respectively. The solitary waves 
established in [14^ [T2] are of elevation and decays exponentially to a horizontal laminar flow 
far up- and downstream. The study of solitary waves of large-amplitude with an arbitrary 
distribution of vorticity remains open. So this arises the question of particle paths in a 
rotational small-amplitude solitary wave. Following the pattern in [6] for irrotational case, 
we prove in this work the corresponding results on particle paths in rotational solitary waves 
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by using the properties of solutions obtained in |14j . We consider in this work only solitary 
waves with negative vorticity, however with modifications our arguments can be applied also 
for positive vorticity. 

Precisely we show that if the vorticity is negative and increasing from bottom to 
surface of the flow, and if the underlying current is favorable, i.e., is moving throughout the 
fluid in the same direction as the wave, then as time goes on the particle moves similarly 
as in the irrotational case [OJ. We also show that if the underlying current is not favorable, 
then in solitary waves with sufficiently small amplitude, some particles above the flat bed 
move to the opposite direction of wave propagation along a path with a single loop or a 
single hump. Note that this single-loop kind of path does not exist in the irrotational case. 
In [13], as in most of the works on waves with vorticity, the author considered only waves 
that are not near breaking or stagnation, i.e., the speed of an individual fluid particle is far 
less than that of the wave itself throughout the fluid domain. We consider only such waves 
as well. We do not consider the case of wave with stagnation which is however studied in 
[6] in the irrotational case. 

The paper is organized as follows. In Section [2] we present the mathematical formu- 
lation for the solitary waves and recall from [14] some useful properties of the established 
solitary wave solutions. Section [3] contains some conclusions on the vertical and horizontal 
velocity that are relevant for our purposes. The main result is presented and proved in 
Section HI In the final section we give two examples which lie in our settings. 



2. Preliminaries 

We first describe the governing equations for rotational solitary water waves and then 
recall some properties available on their solitary wave solutions. 

2.1. The governing equations for rotational solitary water waves. Choose Cartesian 
coordinates (X, Y) so that the horizontal X-axis is in the direction of wave propagation and 
the y-axis points upwards. Consider steady waves traveling at constant speed c > 0. In 
the frame of reference moving with the wave, which is equivalent to the change of variables 
{X — ct, Y) I— )• {x, y), we use 

17^ = {(x, y) : -d < y < r]{x)}, < d < oo, 

to denote the stationary fluid domain and (u{x,y),v{x,y)) to denote denote the velocity 
field, and define the stream function il^ix^y) by ilj{x,rj{x)) = and 

ipy = u- c, = -V. (1) 

Consider also only waves that are not near breaking or stagnation, so that tpy{x, y) < —6 < 
in Clfj for some 5 > 0, which implies that the vorticity u; = Vx — Uy is globally a function 
of the stream function ^jJ, denoted by jii^)', see [8]. The solitary- wave problem is then, for 
given pq < and 7 £ ^^^([0, — po] j I^)) to find a real parameter A, a domain fi^ and a function 
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■0 e C^i^^r,) such that 

V'y<0, {x,y)Gn^, (2) 

AV' = -7W, (3) 

|V^|2 + 2g2/ = A, y = r]{x), (4) 

^ = 0, y = (5) 

= -Po, y = -d, (6) 

7/(2;) —)• as|3;|— )-oo, (7) 

tpxixjy) — 7- as |x| — 7- 00 uniformly for y. (8) 
Here 17 > is the gravitational constant of acceleration, 

_^ i^y{x,y) dy <0 

is the relative mass flux (independent of x) , and the boundary conditions ([7]) and ([8]) express 
that the wave profile approaches a constant level of depth and the flow is almost horizontal 
in the far field, respectively. Moreover we require that the nontrivial solitary wave is of 
positive elevation, i.e., r]{x) > for all x G M. It is therefore symmetric about its single crest 
and admits a strictly monotone wave profile on either side of this crest (see |15]). Assuming 
the wave crest is located at x = 0, the solitary-wave problem ©-([S]) is thus supplemented 
with the symmetry and monotonicity conditions 

ip{—x,y) = i/j{x,y), r]{—x) = rj{x), and r/'(x) < for x > 0. (9) 

We refer to [HUH] for more details on the derivation of the system ©-([S]). 

2.2. Rotational solitary water waves. We collect some properties of the solitary wave 
solutions established in |T3]. Let 

T{p)= \^ -i{-s)ds andrmin= min r(p) < 0. (10) 

Jo pe[po,o] 

Given j?o < and 7 G C^{[Q,—pq\;M.), for each A G (— 2rmin,oo) the system (H])-® admits 
a trivial solution pair (r/(x), ^(y)) defined on Qq, where r]{x) = 0, the stream function ^{y) 
is x-independent and is the inverse of the function 

y{^)=r^ -d, 
Jpo + 2T{p) 

and the fluid domain 

"0 dp 
'P- ^/\ + 2T{p)' 

The corresponding relative horizontal velocity and vertical velocity are thus given by 

U{y)-c = ^y{y) = -^/Y+W{^W)). n^,y) = -^.(y) = 0. (11) 
Note that 'I'(O) = and ^'(-(i) = -po- Thus 

C/(0) = c-^/\ and U{-d) = c - ^J\ + 2T{po). 

Throughout this paper, we adopt the terminology in [[9j to say this underlying trivial flow 
is favorable if U{y) > for all y G [— (i, 0], is adverse if U{y) < for all y G [— d, 0], and 
is mixed if U{y) changes sign. Note that favorable flow moves in the same direction as the 
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wave propagation (i.e., to the right), while adverse flow moves in the opposite direction of 
the wave propagation. 

To ensure the existence of nontrivial small amplitude solitary waves, the parameter A 
must be chosen to satisfy A > Ac but close to Ac, where Ac > — 2rinin is the unique solution 
of 

r° = 1 (12) 

(Ac + 2r(p))=^/2 g 

For each such a given A and for given pQ < 0, it was shown in [H] that there exists a 
nontrivial small amplitude solitary-wave solution pair (ri{x),tlj{x,y)) to ©-([H]) defined on 
Clyj, with r]{x) satisfying 

\rj{x)\ + \rj'{x) \ + |r?"(x)| < (A - Ac)r for all x G M (13) 

and for some constant r > independent of A, and the corresponding horizontal velocity 
satisfying the following properties: 

(iE^t), ri{x)) — )■ c — \/A = f/(0), as |x| — )• oo; 
(zE(2;),7/) —7- U{y) as — )• oo uniformly for y. 

The first property is obvious since by ([7]) and dS]), we have {u{x, ry(x))— c)^ = ^pyi^x, r]{x)) — ) 
A as \x\ — oo, which equivalently gives (F[T]) as we have assumed u{x, y) — c = ipyi^i y) < in 
O^, while the property (F|2j) can be deduced from the construction of solitary wave solutions; 
see [Hj. Indeed for A = Ac + e with e > 0, there exists a function w^{q,p), whose derivatives 
with respect to {q,p) up to order 2 tend to uniformly for p as \q\ — )• oo, such that the 
horizontal velocity is determined by 

u(^x y^ = c — 

(A + 2r(-V'(x, y))-'/' + ewl,{V^x, -i^{x, y)) ' 

where Wp{q,p) denotes differentiation in the p- variable. 

We conclude this section by recalling some properties of streamlines. Due to 'i/'y < 
throughout (iri, we have that for all p G [— P0)0] the streamline 

{{x,y) : ij{x,y) = -p} 

is a smooth curve y = crp{x). Note that 

ao{x) = ri{x), ap„[x) = -d and cr (x) = — — — -. (14) 

^ '4jy{x,ap[x)) 

Observing the fact that = is bounded and that ipx{x-,y) — >■ uniformly for y 
as |x| — )• oo, we deduce that for each fixed yo G [Oi'^(O)] the streamline y = (Tp(x) with 
p = — -0(0,7/0)1 passing through the point (0, yo)) has an asymptote y = l{yo) as |x| — )• 00, 
with /(r?(0)) = and l{-d) = -d. 

3. Vertical and horizontal velocity 

We first divide the fluid domain $7^ into two components 
Sl_ = {(x, y) G : X < 0, < y < r]{x)} and = {(x, y) G : x > 0, < y < r/(x)}, 
and denote their boundaries by 

S'_ = {(x, y) G : X < 0, y = r/(x)}, 5_ = {(x, y) G : x < 0, y = -d}, 
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respectively 

S+ = {{x,y) e M2 : X > 0,y = B+ = {{x,y) G : x > 0, y = -d}. 

Lemma 3.1. Suppose that 7'(p) < for all p G [0, \po\]. Then 

(a) v{x, —d) = for all x € M, and t;(0, y) = for y G [— d, /^(O)] ; 
f6j f (x, y) < if (x, y) G U S"-, and i;(x, y) > if (x, y) G $1+ U 5+; 
(^cj Vy{x,—d) < i/x < 0, and Vy{x, —d) > i/x > 0; 
f,(0,y) >0/ory G (-(i,r/(0)). 

Proof. Since = — V'x) the first result follows from ([6]) and the symmetry property ^{—x, y) = 
tpix,y) in ©. 

Next we only prove the lemma for x > and the results for x < can be proved 
similarly. Differentiating ([5]) with respect to x gives v = —^x = from which follows 

v > for (x,y) £ S+ as ipy < and 7/'(x) < for x > in view of ([2]) and ([9]). To prove 
> in we assume first on the contrary that there exists a point (xo,2/o) £ ^+ such 
that v(xo,yo) = — e < 0. Then we can find a bounded domain 

= {(x, y) G M : < X < A;, -d < y < r/(x)}, k G M+, 

such that (xo,yo) G ^^+,fc- Moreover in view of ([8]) we can choose k sufficiently large so that 
v{k,y) > —e for y G [—d,r]{k)]. This means that v attains its minimum at the interior point 
{xo,yo) of which contradicts to the strong maximum principle applied to v on the 

domain J^+^fc, as v satisfies Av + ^'{'ip)v = with 'y'{p) < by differentiating ([3|). Therefore 
we have v > in If v = at a point (xo,yo) of Again we can choose a bounded 
domain il^ ^ containing (xo,yo)- Then t; > on the boundary of 0^ ^. We can thus apply 
the strong maximum principle on ^ once more to conclude that u = on ^ , which 
contradicts w > on the half surface 5"+. This proves f; > in 

Since v attains its minimum in Ctj^ on the half bed Bj^ and on the crest line x = 0, 
Hopf's maximum principle ensures that > on Bj^ and Vx> ^ on {(0, y); —d < y < ??(0)}, 
completing the proof. □ 

The above lemma combined with (jl4|) and the fact that ipy < gives 

Corollary 3.2. The streamline y = ap{x) with p G (po;0] satisfies that cr'p{x) > for x < 
and o"p(x) < for x > 0. 

Lemma 3.3. Suppose that j{p),j'{p) < for all p G [0, \po\], that < \ — Xc < e with e 
small such that nontrivial solitary wave exists, and that c > + 2T{pq). Then u{x, y) > 
for X £ Q,n- 

Proof. We assume throughout this proof that 7(p) ^ for p G [0, \po\], since we already know 
that for 7(p) = 0, i.e., the irrotational case, u{x,y) > in il^; see [6]. Recall that U{y) is 
the horizontal velocity of the trivial laminar flow. Thus U'{y) = ^yy{x,y) = — 7(^'(y)) > 0, 
and there exists an interval I such that the strict inequality holds for y £ I since 'y{p) ^ 
by assumption. This combined with and the assumption c > + 2T{pq) gives 

U{y) > U{-d) = c-yj\ + 2T{pq) > for y G [-d,0], and C/(0) > 0. 
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It was proved in [T7] that if 7(p) < for all p E [0, |po|]) then 

—u(x,ri{x)) > for X < 0, and —u(x,r](x)) < for x > 0. (15) 
ax ax 

In other words, along the free surface u increases from x = — oo to the crest x = 0, and 

thereafter it is decreasing. On the bottom S = {(x, y) G : y = —d}, we have, in view of 

Ux = —Vy = ipxy and Lemma [XT] - fcl). 

Ux{x, —d) > for X < 0, and Ux{x, —d) < for x > 0. (16) 

Then the fact that U{—d) > and U{0) > together with the monotonicity properties (|15p 
and ()16p yields u > on the free surface and on the bottom. Differentiating ([3]) with respect 
to y yields that u satisfies 

Au + j'{tp)u = 7'(V')c < 0. 

Finally considering U{y) > for y G [— d, 0] and the properties (F[T])-(F|2]), we can deduce 
tt > in fi^ by using the strong maximum principle as in the proof of Lemma l3.H -(jb]). 
completing the proof. □ 

The above lemma considers the case c > + 2r(po)i that is the underlying flow is 
favorable. For the remaining cases, we can only derive some conclusions on the horizontal ve- 
locity u for some special vorticity function 7 and for A sufficiently close to the corresponding 
Ac defined in (jl2p . To be exact, we prove 

Lemma 3.4. Let Ac be determined in U2\) by the vorticity function "f{p) with p G [0, \po\]- 

(a) // 7'(p), 7"(p) < for p G [0, \po\] and 7(0)-v/A^ > —g, then there exists ei > such that 
for A G (Ac, Ac + El), we have 

Ux > for (x,y) G and < for (x,y) G ^1+; 

(b) If j{p) < for all p G [0, \po\], then there exists £2 > such that for A G (Ac, Ac + £2); 
we have 

Uy > for (x, y) G Cl. 

Proof, daj) Similar results have been obtained for periodic waves in [10], so we will adapt 
the proof there to our present solitary wave case. We prove only the result for (x, y) G 
and the rest can be proved similarly. Note it suffices to prove that the stream function ■0 
satisfies 

Ipxy < for (x,y) G 

Since ip{x,ri{x)) = 0, we have ipx = —il'^^y the surface, which combined with ^ 

gives 



Mx,r]{x)) = -V(A-2OT)/(l+r/'2). 
Differentiating the above equality and ipx = ~'n''^y along the surface gives respectively 



i^xy + v'i^yy = -dx\ . ,2 and ipxx + v'i'yy + 2??Vxy + ??'V?/ = 0, on y = rj{x). 



'X-2grj 
1 + r?' 

Moreover on the surface one has 

Ipxx + Ipyy = -7(0) 
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Combination of the above three equahties gives 

^'(5(l_^'4) + 2V'(A-2ffr?)+7(0)^/A^2^(l + V2)3/2) 

= VA - 2g,(l + • 

In view of ()13p and 7(0) -v/A^ > —5; we have when A is sLifficiently close to Ac that 



g(l _ n'^) + 2ri"{\ - 2ot) + 7(0)VA-2ot(1 + t?'^)^/^ > 0. 

This gives ipxy < on the surface since ry' < for x > and — 2gr]{l + r/'^)^/^ > 0. Since 
v{0, y) = 0, we have ipxy = —Vy = on the hne x = 0. On the bottom ilixy < holds due to 
Lemma Ism ffcj). And tpxy = — )• as x — )• 00 since u{x,y) —)•[/(?/) as x 00. Finally ^xy 
satisfies 

(A + j')i;xy = -l>xi>y > 0. 

The conclusion follows from the maximum principle. 

Recall from [14.J that for A = Ac + e, there exists a function w^{q,p) such that 
u{x, y) is determined by 

" " (A + 2r(-^(x, y)))-i/2 + ewl,{^ex, -^{x, y)) ' ^^'^^ 

where Wp{q,p) stands for the differentiation of the function w^{q,p) with respect to p, and 
the family {w^{q,p);0 < e < 1} satisfies that for e small 

sup{\dp^w'{q,p)\-j + k<2,{q,p)eRx [po,0]] < r, (18) 

with r > some constant independent of e. Differentiating (I17p with respect to y gives 

-i^y {-iW (A + 2r(-V^))-3/2 + ew^ppi^x, -^)^ 



((A + 2r(-V))~^/' + ew^^x, -^))' 



Since 7(p) < 0, if we denote 7max = maXpg[ojp,,|] 7(p), then 7max < and Tmax = maXpg[p^^ q] > 
0. As a result, observing 

-7(^) (A + 2r(-^))-3/2 + eu;;p(\/^x, -V) > -7max (A + 2rmax)"^/^ + ewl^i^/Ix, -t^) 
and 

lim (-7max (A + 2rmax)"^/^ + ew'ppiV^X, -V-)) = -7max (Ac + 2rmax)"^/^ 

due to (fT8]l . we get, in view of — 7max (Ac + 2rinax)^^^^ > 0, that 

-7('0) (A + 2r(-V'))"^/^ + ew'ppiVIx, -4^) > 

when £ is sufficiently small. This combined with —ipyyO gives Uy > for e sufficiently 
small. □ 

Lemma 3.5. If ^{p) < 0,j'{p),j"{p) < for p G [0, |po|] o,nd 7(0)\/A^ > —g, then for 
A G (Ac, Ac + £0) with Eq = min{ei,e2}! along every streamline y = ap{x) with p £ {pQ,0), 
the horizontal velocity u is strictly decreasing in and strictly increasing in as a 
function of x. 
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Proof. Since 

(T„(x)) = Ux + Uya Jx) = Ux + Uy—^ 
dx y ' " u — c 

due to (|14p , the conclusion follows immediately from Lemma 13.11 and Lemma 13.41 □ 

4. Main result 

Recalling that y = l{yo) is the streamline asymptote introduced at the end of Section 
[21 we are now ready to state our main result describing the particle trajectories in a solitary 
wave with negative vorticity. 

Theorem 4.1. Let the flux po < 0, the vorticity function 7 € C^(0, |po|) be given, and 
let Ac be determined as in [W\) . Assume that j{p) < 0,7'(p) < for all p G [0, |po|]; o.iT'd 
that A > Ac such that nontrivial solitary wave solutions to ^-^Si) exist. Then the following 
results hold. 

(a) Any particle above the bed reaches at some instant to the location {Xq, Yq) below the wave 
crest (Xo, T?(0)); 

(b) For c > ^X + 2T{po), as time t runs on (—00, to); the particle above the flat bed moves 
to the right and upwards, while for t > to the particle moves to the right and downwards, 
as in Figure (a); the particle on the fiat bed moves in a straight line to the right at a 
positive speed; 

(c) If assume additionally that 7(0)-v/A^ > —g o-nd ^"{p) < for all p £ [0, |po|]; then there 
exists eo > such that for A € (Ac, Ac + eo) o,nd 

(i) for y/X < c < ^A + 2T(j)q), there does not exist a single pattern for all the particles 
above the flat bed; ifu{x, —d) > 0, depending on the relation between the asymptote 
y = 1{Yq) and the zero point y^ ofU{y), some particles move to the right along a 
path with a single hump as described in some move along a single-loop path to 
the left, as in Figure{l\-(b); if u{x, —d) < 0, there are additionally some particles 
moving to the left along a single-hump path; see FigureU\-(c); 

(a) for c < y/X, depending on the signs o/u(0,r/(0)) and u{0,—d), there are three 
possibilities for the particles above the flat bed; see Figure 0' 

(Hi) for a particle on the flat bed in both the cases (i) and (ii), if u{x, —d) < it moves 
to the left in a straight line, while if u{x, —d) > it firstly has a backward-forward 
pattern of motion and then moves to the left in a straight line; 

(d) The particle trajectory is strictly above the asymptote Y = 1{Yq) of the streamline Y = 
ap{X — ct) with p = — ^'(O, Yq)- 



Proof. The path (past and future) {X{t),Y{t)) of a particle with location {X{0),Y(0)) at 
time t = is given by the solution of the non-autonomous system 

r X = u{X -ct,Y), 
\ Y = v{X - ct,Y). 

Working in the moving frame x = X — ct and y = Y, we transform the above system into 

( x = uix,y) - c, 
\ y = v{x,y). 
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u(0, -d)>0 u(0, -d)<0 



Figure 2. Particle path above the flat bed in a smaU solitary wave with a mixed un- 
derlying current. 

This is a Hamiltonian system with Hamiltonian 'tp{x,y) in view of ([1]), meaning that the 
solutions {x{t),y{t)) of ([IHD he on the streamhnes. Ah solutions of (fTn|) are defined globally 
in time since the boundedness of the right-hand side prevents blow-up in finite time. 

(jaj) Since in the moving frame the wave crest is assumed to be located at x = in our 
setting, the sign of x{t) describes the position of the particle with respect to the wave crest 
at time t: the particle is exactly below the crest if x(t) = 0, is ahead of the crest if x{t) > 
while is behind the crest if x{t) < 0. Note that x = u{x,y) — c < —5 < throughout fl^. 
This uniform upper bound on x implies that x{t) — )• =poo as t — t- ±00 and there is a unique 
time to such that x(to) = 0. That is, to each fiuid particle moving within the water there 
corresponds a unique time to ^ so that at t = the particle is exactly below the wave 
crest, while afterwards it is located behind the wave crest, the wave crest being behind the 
particle for t < to- 

In the subsequent of the proof, we always assume that the particle is located below the 
crest at time t = to at the location {Xq,Yq) = (X(to), ^(to)), or equivalently (j;(to), y(to)) = 
(0,10) which implies in fact Xq = cto- 

(|b]) Assume c > -y/A + 2T{pq). For a particle located above the fiat bed, we prove 
only the statement for t > to and the case t < to can be proved similarly. Since x(t) is 
strictly decreasing, one has X{t) — ct = x{t) < for t > to, which combined with Lemma 
l3.H -(jb|) implies y < for t > tQ. Furthermore we have X > for all time by Lemma 13.31 
That means the particle moves to the right and downwards as time runs on (to, +00). The 
statement for particles on the bed follows from Lemma lS.lK faj) and Lemma [3. 31 The particle 
path for this case is depicted in Figured]- (a). 

(jci)-(i) Since VX < c < ^JX + 2r(po), we have C/(0) > and U{-d) < 0, thus U{y) has 
a unique zero point, say at y*. Assume first that m(0, —d) > 0. The monotonicity properties 
of u provided by Lemma l3.4l -f[ai) and ()16p show that in Cl^ the level set {u = 0} consists 
of a continuous curve C+ in the moving frame. The curve is confined between y = —d and 
y = y^f, and can be parameterized by x = h{y) with h'{y) > 0, h{y) — )• -|-oo as y ^ y^, and 
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u(0, -d))>0 «(0, -d) < < 11(0,77(0)) 11(0, -d) < 0,w(0,r;(0)) < 



Figure 3. Particle path above the flat bed in a small solitary wave with an adverse 
underlying current. 

h{—d) > with the equality holding when u{x, —d) = 0. In r2_ the level set {u = 0} is given 
by the reflection C_ of the curve C+ across the line x = 0. Below the curve C+ and C_ we 
have u < (including the bottom), while above and between the two curves we have u > 0. 
Furthermore, in view of Lemma 13.5^ if a streamline and the curve Cj^ (rep. C_) intersect, 
they intersect exactly once. 

Recall that {x{t),y{t)) lies on the streamline y = cjp(x) with p = —Tp{0,Yo) since 
{x{to),y{to)) = (OjYq), and recall that y = 1{Yq) is the asymptote of the streamline passing 
through the point (0, Yb)- In virtue of Corollarv 13.21 the path {x{t),y{t)) is located below 
y = Yq and above the asymptote y = 1{Yq) for all time t. If 1{Yq) > y^, then we have from 
the above arguments that u > for all time t, while < if t > to and f > if t < to- 
This is the same situation as described in (jb]), so that the particle moves to the right along 
a path with a single hump, as depicted in Figure [!]■( a). 

If 1{Yq) < 2/^,, the particle trajectories above the flat bed are as shown in Figure [T]-(b). 
Indeed, in this case, as time t increases from —00, the path {x{t),y{t)) intersects successively 
the curve C+ from below at t = t+, the vertical line x = from right at t = to, and the 
curve C- from above at t = t_. In the time interval t G (—00, t+) U (t_, +00) we know that 
u < so that in the physical variables {X,Y) the particle {X{t),Y{t)) moves to the left. In 
the time interval t £ (t-(_,t_) we have n > so that {X{t),Y{t)) moves to the right. Also 
we have v > when t < to so that {X{t),Y{t)) moves up, while when t > to we have v < 
so that (X(t),y(t)) moves down. Thus in this case the particle above the flat bed moves to 
the left along a path with a single loop. 

It remains to treat the case u{0, —d) < 0. Note ti(0,r/(0)) > in virtue of U{0) > 
and fflS]) . Thus, in view of Lemma [3^(jb]l . there exists a unique y £ {—d,r]{0)) such that 
u{0,y) = 0. Observe y < y^: since ti^: < in r2+ and u{x,y) — C/(y) as x — )• 00. The 
corresponding curve C+ is now located between y = y* and y = y, intersecting the line x = 
at y = y. For Yq > y, the paths {X{t),Y(t)) are similar as encountered in the case when 
u{0, —d) > 0. While for Yq < y, we have {x{t),y{t)) is below C+ and C_, so that n < for 
all the time. As a result, the particles move to the left along a single-hump path; see Figure 

mic). 

We depict all the possible trajectories in this case in Figure [2j 

(jci)-(ii) This case can be treated similarly as that in the above, so we omit the details 
and show the particle trajectories in Figure [31 

(jc])-(iii) The conclusions for particles on the flat bed follow from Lemma |3. 11 ffaj) and 

m- 
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(pl Observing Y = v{X — ct,Y), X{t) — ct = x{t) — )• — oo as t — )• +00 and v{x,y) 
converges to as — )■ 00 uniformly in y, we have the existence of some a such that 
hmj_>.+oo ^(i) = Oi. Since ip{x,y) is the Hamiltonian function of the system (jl9p . we have 

iP{x{t),y{t)) = ^P{x{to),y{to)) = i^iO,Yo). 

Thus y{t) = (jp(x(t)) with p = —■iP{0,Yq). Consequently 

lim Y{t) = lim y{t) = lim crp{x{t)) = l{Yo). 

The proof is completed. □ 

Remark For the path with a single loop, if we define the size of the loop by the net 
horizontal distance moved by the particle between the two instants t+ and t_ when its 
horizontal velocity changes sign, that is, 

X{t-)-Xit+), 

then the size decreases with depth. In fact it can be computed by 

ri- dX ft- 
X{t^) - X{t+) = \ —dt=\ u{x,ap{x))dt. 

Since Uy > and ^^^^^ = — > 0, we have X{t-) — X{t^) decreases as p decreases. 

5. Examples 

In this final section we give two examples of the vorticity function which satisfy the 
conditions imposed in our main theorem. 

5.1. Negative constant vorticity. In the case of negative constant vorticity ^{p) = — 
for p G [0, |po|] with ujq > 0, we have obviously 'y'{p),j"{p) < 0. So it remains to verify 
7(0)\/A^ > —g, or equivalently, Ac < (/^/wq. Recall from ([12]) that Ac > — 2rinin is such that 

rO dp _ 1 

■^P" (Ac + 2r(p))^/2 ~ 9 
holds, where T{p) = —ujqp and Fmin = in this case by (fTOl) . Set 

-0 dp 



F(A) = r 

J Pi 



'p« (A - 2ujopf^^ ' 

Then direct computation shows that F'{X) < 0, F{X) — )• +00 as A — )■ 0+, and 

1 1 1 

hm F(A = , < -. 

A->3V<^j5 9 V^^^wgpo 9 

Thus there exists a unique Ac G {0,g'^/iOQ) such that F{Xc) = l/g. 

5.2. Negative afRne linear vorticity. For the vorticity j{p) = —ap + b with a > and 
6 < 0, we have obviously ^{p) < 0, 7'(p),7"(p) < for p € [0, \po\] and 

= r — = r — — ^, 

Jp« (A + 2F(p))^/^ ■^P" (A + ap2 + 2bpf'^ 

with A > — 2Fmin = 0. Then we may verify as above that -F(A) is strictly decreasing with 
limA-j.0 -^(A) = +00 and Iim.x^g2/h2 F{X) < 0. This gives the existence of Ac G (0,5(^/6^) 
such that -F(Ac) = 1/5 and consequently 7(0) -v/A^ > —g- 
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